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Quantum creation processes during the very rapid early expansion of the universe are 
believed to give rise to temperature anisotropies and polarization patterns in the CMB 
radiation. These have been observed by satellites such as COBE, WMAP, and PLANCK, 
and by bolometric instruments placed near the South Pole by the BICEP collaborations. The 
expected temperature anisotropies are well-confirmed. The B-mode polarization patterns in 
the CMB are currently under measurement jointly by the PLANCK and BICEP groups to 
determine the extent to which the B-modes can be attributed to gravitational waves from 
the creation of gravitons in the earliest universe. 

As the original discoverer of the quantum phenomenon of particle creation from vacuum 
by the expansion of the universe, I will explain how the discovery came about and how 
it relates to the current observations. The first system that I considered when I started 
my Ph.D. thesis in 1962 was the quantized minimally-coupled scalar field in an expanding 
FLRW (Friedmann, Lemaitre, Robertson, Walker) universe having a general continuous scale 
factor aft ) with continuous time derivatives. I also considered quantized fermion fields of 
spin-1/2 and the spin-1 massless photon field, as well as the quantized conformally-invariant 
field equations of arbitrary integer and half-integer spins that had been written down in the 
classical context for general gravitational metrics by Penrose. 

It was during 1962 that I proved that quanta of the minimally-coupled scalar field were 
created by the general expanding FLRW universe. This was relevant also to the creation of 
quantized perturbations of the gravitational field, since these perturbations satisfied linear 
field equations that could be quantized in the same way as the minimally-coupled scalar field 
equation. In fact, in 1946, E.M. Lifshitz had considered the classical Einstein gravitational 
field in FLRW expanding universes and had shown that the classical linearized Einstein held 
equations reduced, in what is now known as the Lifshitz gauge, to two separate classical 
minimally-coupled massless scalar held equations. These held equations of Lifshitz, when 
quantized, correspond to the held equations for massless gravitons, one equation for each of 
the two independent polarization components of the spin-2 massless graviton. I will discuss 
this further in this article 3 

1 Plenary Lecture given September 2, 2014 at the ERE2014 Conference in Valencia, Spain. 
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I. INTRODUCTION 


In 1962, at Harvard when I began my Ph.D. thesis, I wanted to work at the interface of general 
relativity and quantum field theory. I had already studied general relativity quite carefully as an 
undergraduate before coming to Harvard. I had the good fortune to learn quantum field theory and 
particle physics at Harvard from Wendell Furry, Roy Glauber (Nobel Prize 2005), Sidney Coleman, 
Sheldon Glashow (Nobel Prize 1979), and Julian Schwinger (Nobel Prize 1965). 

I wanted to find new consequences of the quantum field theory of elementary particles in the 
context of Einstein’s theory of general relativity. At the time, I felt that quantizing the nonlinear 
gravitational field itself was so difficult that I would not be able to make significant progress in 
trying to go beyond the deep work that had already been done in that area. Nevertheless, I felt 
that it would be valuable to study quantized elementary particle fields in the curved space-times 
that were solutions of the classical nonlinear Einstein gravitational field equations. Luckily, Sidney 
Coleman agreed to be my thesis advisor on such a project, which was outside the main stream of 
the time. 

I started by looking for new consequences of quantum field theory in the isotropically expanding 
cosmological space-times that were solutions of the nonlinear classical equations of general relativ¬ 
ity. To investigate the consequences of quantum field theory in the expanding universes of general 
relativity, it was first necessary to extend the known quantum field theory of elementary particles 
from the well-established flat Minkowski space-time of special relativity to the context of a classi¬ 
cal general relativistic isotropically expanding universe with a general expansion scale-factor, a(t). 
Some examples that I studied were the dust-filled universe, the radiation-dominated universe and 
the exponentially-expanding universe of deSitter, but my main result was for general a(t). 

The quanta (particles, scalar perturbations, and tensor gravitational wave perturbations) that 
reach a steady-state density during the inflationary stage of the expansion are created or maintained 
by the quantum field theory process of particle creation that I discovered and will discuss in this 
lecture. The quantum field theory vacuum state in the exponentially expanding universe is known 
as the Bunch-Davies[l| vacuum. That vacuum state is symmetric under the deSitter invariance 
group of the exponentially expanding FLRW universe. The Bunch-Davies vacuum is the natural 
vacuum state having the 10 deSitter group symmetries',^], just as the natural vacuum state in 
flat Minkowski spacetime has the 10 Poincare group symmetries. 

The constant, non-zero, average number density of particles, scalar perturbations and/or tensor 
gravitational wave perturbations is consistent with the result of spontaneous particle creation from 
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vacuum of an FLRW universe that undergoes the deSitter exponential expansion. I will discuss this 
further in Section uni and will show that the creation rate of these perturbations and/or particles 
follows immediately from the result that I obtained in my Ph.D. thesis for general continuous 
expansion factors aft) having a sufficient number of continuous time derivatives. This follows very 
simply from my general equation for the average number of particles created, and from the deSitter 
invariance of the exponentially expanding universe. 

I chose to first consider the spatially-ffat FLRW universes because I felt that one could dy¬ 
namically propagate the experimentally verified theory of quantized fields from flat spacetime into 
the curved spacetime of the spatially-flat isotropically expanding universes. It was not necessary 
in this framework to go beyond the already tested assumptions of classical general relativity and 
of quantum field theory in flat spacetime. The dynamical equations of the quantized fields were 
taken to be the simplest generally-covariant ones obtained from the flat-spacetime quantized field 
equations. It was then possible to use the field equations in the curved spacetime to propagate the 
fields and their commutation relations into curved spacetime. Any fundamentally new results com¬ 
ing from this minimal set of assumptions would have to be taken very seriously. As it turned out, 
this procedure also gave an independent new proof of the spin-statistics theorem; for bosons only 
the commutation relations were propagated consistently from the initial to the final flat spacetime, 
and for fermions only the anticommutation relations were propagated consistently. 

First, I considered the quantized generally-covariant minimally-coupled free scalar-field equation 
in an isotropically-expanding spatially-flat universe with a continuous and smooth scale factor aft). 
Because there was no ambiguity as to the definition of the measurable particle number operator 
of the quantized field in flat spacetime, it was natural to consider what happened to the particle 
number when the scale factor aft) went smoothly from an initial constant value a\ at early times 
to a final constant value 02 at late times. During the periods of time when a± and 02 were 
constant, simple constant rescalings of the spatial coordinates at early and late times showed that 
spacetimes before and after the expansion of the universe were Minkowskian, each exhibiting the 
10-dimensional Poincare group of symmetries. The number density of particles of the quantized 
field was taken to be 0 at early times. Evaluating the expectation value of the number density at 
late times, after the universe had undergone a smooth expansion, I found that the number density 
of created particles was finite and nonzero! This demonstrated that real particles were created 
from the initial vacuum state by the expansion of the universe. 

To briefly summarize the calculation, I used the Heisenberg picture to evolve the quantized field 
during the expansion of the universe. The state vector of the quantized field was chosen to be the 
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well-defined vacuum state in the initial Minkowski spacetime. As a result of the time-dependent 
unitary evolution of the quantized field, the creation and annihilation operators of the quantized 
field at early times had evolved into linear combinations of the creation and annihilation operators 
of the quantized field at late times. Consequently, the state vector having no particles of the 
quantized field at early times does have a non-zero density of particles of the quantized field at 
late times. This implies that particles are created by the expansion of the universe. 

My results first appeared in my Ph.D. Thesis 4] (“The Creation of Particles by the Expansion 
of the Universe”, Harvard University, 1966) and also in [f]-7]. My thesis also included results I had 
obtained for the spin-1/2 fermion field, including particle creation for fermions of non-zero mass, 
as well as the absence of particle creation for the massless Pauli neutrino field. Those results for 
fermion fields appeared in [7]. 

I also found that particles of 0 -mass satisfying the minimally-coupled scalar field equation would 
be created from the initial vacuum state as a result of the expansion from an initial Minkowski 
spacetime to a final Minkowski spacetime. This result implied that quantum gravitons would be 
created by the quantized linearized Einstein gravitational held equations. E.M. Lifshitz 8] in 1946 
had published an important paper on the classical linearized Einstein gravitational held equations, 
in which he proved that there is a gauge (the Lifshitz gauge) in which the linearized Einstein held 
equations reduce to a set of two separate minimally-coupled scalar held equations. One could then 
use the Lifshitz gauge to directly apply my results to the Einstein gravitational held to show that 
quantized gravitons (tensor perturbations) were created by the expanding universe. A paper on 
this was published by Leonid P. Grischukp], in which he used the Lifshitz gauge to apply my earlier 
quantization of a minimally-coupled scalar-held, including my expression for particle creation in a 
general expanding universe. He used the Lifshitz gauge to obtain an expression for the creation of 
gravitons in a universe with a{t ) proportional to a power of t. 

A more detailed quantization of the linearized Einstein gravitational held was published in 
1977 by Lawrence H. Ford and me in [bjj]. One of the things we showed in our paper was that 
there were spurious infrared (IR) divergences present in Grishchuk’s result that were not present 
when we imposed more appropriate late-time boundary conditions. In our paper, we also reviewed 
the derivation of the Lifshitz gauge and performed an explicit quantization, defining creation and 
annihilation operators for massless gravitons in each of their two polarization states. (We carried 
this out for more than one way of specifying the pair of gravitational wave polarization states.) 

The minimally-coupled scalar held equations that appear in the Lifshitz gauge are the same 
ones I studied 4j-[6j in the 1960’s for general a(t). In the exponentially expanding inhationary 
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universe, there are well-known analytic solutions of these equations that can be applied to the 
creation of gravitons and gravitational wave perturbations present in the inflationary universe. 
These form the basis for obtaining the inflationary universe predictions of linearized gravitational 
waves and gravitons created in the early universe, as well for scalar gravitational perturbations 
that are related to the temperature variations of the CMB. 

In Section mi will show how one obtains directly from my published results of the 1960’s, the 
average number of quantized particles created as a function of time in an expanding physical volume 
of the exponentially expanding universe. This result follows directly from deSitter invariance of 
the exponentially expanding universe. In addition, if there are similar particles already present 
at the time that the exponential expansion begins, then the initial presence of these particles can 
significantly increase the rate of particle creation. This stimulated creation of identical bosonic 
particles occurs also for general a(t). (For fermions the initial presence of identical fermions would 
decrease the rate of particle creation.) 

Scalar perturbations of the gravitational field were also created by the exponential expansion 
of the universe, before reheating occurred. Reheating resulted in the formation of a plasma of 
charged particles in a high- temperature, radiation-dominated expanding universe. The acoustic 
waves that these scalar gravitational perturbations produced in the plasma were present for about 
400,000 years, until the universe cooled enough for electrons and protons to form neutral atoms, 
after which the photons of the CMB were able to propagate freely (for the most part). However, 
temperature differences in the plasma acoustic waves present at the time of last scattering resulted 
in temperature differences in the free-streaming CMB radiation. These temperature differences in 
the CMB radiation at the time of last scattering were responsible for the temperature anisotropies 
of about 1 part in 10 5 that were ultimately measured by the COBE satellite. 

If the recent possible gravitational wave observations are confirmed and if the intensities of the 
B-mode polarization patterns they produce are consistent with the inflationary universe model, 
including the observed temperature anisotropies in the CMB, then this consistency would be further 
evidence for an early inflationary expansion of the universe. 

In the next Section, I will briefly explain how I proved in my Ph.D. thesis (see also Bfl) 
that minimally-coupled scalar particles are created by expanding FLRW universes. 

But first let me quote the Conclusion of my 1968 paper [5]: “The particle creation in the expand¬ 
ing universe at the present time is quite negligible. However, for the early stages of a Friedmann 
expansion it may well be of great cosmological significance, especially since it seems inescapable if 
one accepts quantum field theory and general relativity. In considering the large amount of par- 
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tide creation taking place in the early stages of an expansion, it is necessary to take into account 
the reaction of the matter created back on the gravitational field. Furthermore, it may be neces¬ 
sary to consider the effects of the quantization of the gravitational field. Therefore, no conclusive 
quantitative result can yet be reported here concerning the primeval creation.” 

This conclusion of my paper emphasized the strong effects that particle creation would have on 
the early expansion of the universe, so the reaction-back would have to be taken into account. The 
first such calculation was done by me and Steven A. Fulling jll], in which we showed by numerical 
evolution of the semi-classical Einstein equations that a bounce could occur as a result of the 
particle creation process, thus avoiding the cosmological singularity. Our result also demonstrated 
that the classical energy-conditions can be avoided through the creation of particles by quantum 
field theory in an FLRW universe. 

. n 

In my thesis |4[, I also used quantum measurement theory to develop the method of adiabatic 
regularization (or renormalization) as applied to the particle number operator during the expansion 
of the universe. In 1974, Fulling and I extended the adiabatic regularization method to renormalize 
the expectation value of the energy-momentum tensor 0,Q,Q. 

An interesting sidelight, is that I was hired as an Instructor to teach at the University of North 
Carolina (UNC) for 1966 and 1967. My thesis advisor was away in Europe, so I had to wait until he 
returned to Harvard in 1966 to fly to Boston for my thesis defense. The thesis defense committee 
consisted of Sidney Coleman, Sheldon Glashow, and Walter Gilbert (Nobel Prize 1980). 

While teaching at UNC, I attended a Colloquium given at UNC in 1966 or 1967 by Fred Hoyle. 
After his talk, I explained to him my as yet unpublished work. He was an advocate of the steady 
state theory (of Bondi, Gold, and Hoyle), which assumed a permanent exponentially expanding 
universe having a constant density of particles, maintained by a process that they had postulated, 
in which particles or atoms were constantly being created, including at the present time, at a 
sufficient rate to maintain a steady-state universe that did not change significantly on cosmological 
time scales. I explained to him that my quantum field theory mechanism of particle creation 
would not support such a steady-state universe at the present time because the particle creation 
rate would be much too small. However, I explained to him that my process of quantum field 
theory particle creation could support an exponentially expanding steady-state universe at very 
early times with a much higher particle density. He listened, but was not enthusiastic about the 
idea because at that time he believed in a universe that was eternal and unchanging. 

In Section 1IIII of this article, I will show how the equations that I first derived in my 
Ph.D.thesisjj-jL] for the general case of particle creation by an FLRW universe directly yields 
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the rate per unit physical volume at which particles or quanta are created in an exponentially ex¬ 
panding, spatially-flat FLRW universe that is in (or has evolved to) a deSitter invariant state. For 
example, the inflationary universe is in such a state prior to reheating. This very high rate of parti¬ 
cle creation is what I was referring to in connection with the exponentially expanding steady-state 
universe in the discussion I had with Fred Hoyle after his Colloquium at UNC. 


II. MINIMALLY-COUPLED SCALAR FIELD IN FLRW UNIVERSE 


Let us now consider the case of a free quantized scalar field <f> in a smoothly changing spatially 
flat FLRW universe with line element 


and equation of motion 


ds 2 = dt 2 — a 2 (t)(dx 2 + dy 2 + dz 2 ), 


(□ + = 0 , 


(1) 


( 2 ) 


where □ = g ,w V tJ V u . with Va denoting the covariant derivative. 

For now we are considering the so-called minimally-coupled scalar field of mass m (which may 
be 0). More generally, one could include a coupling to the Ricci scalar curvature, R = g^R^, 
namely, (□ + m 2 + £R)4> = 0, where £ is a dimensionless constant. We work in units with h and 
c each equal to 1. (For brevity, I will not go into spin-1/2 and spin-1 fields, and spatially-curved 
FLRW universes, which I also considered in my Ph.D. thesis[4-7|. 

Writing ([2]) with the metric of dTJ) , one finds that 


a 3 dt(a 3 dt4 >) — a 2 ^ d 2 (j) + m 2 (j> = 0. 


( 3 ) 


It is convenient to impose periodic boundary conditions in a cube having sides of coordinate length 
L and coordinate volume V = L 3 . The physical volume of this cube is expanding, with the 
physical volume proportional to a(t) 3 . As in Minkowski spacetime, this is a mathematical device, 
with L taken to infinity after physical quantities such as the density of created particles have been 
calculated. Then we can expand the field operator cj) in the form 


where 


t { A kk& *) + *)} > 


ft = {2Va{t) 3 )~ l ' 2 e its h k (t). 


( 4 ) 


( 5 ) 
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Here k l = 2itin l /L with n l an integer, k = \k\, and hf.(t) satisfies 

^ 2 hk(t) + ~ 2 ^k(t) + m 2 h k {t ) - ^ (^j h k (t ) - ^ L h k(t) = 0. (6) 

Examples in which this equation has analytic solutions are power law expansions and the ex¬ 
ponential expansion. For the latter case, a(t) oc exp (Ht) with H being the Hubble constant during 
the exponential inflationary stage of the expansion, the solutions are Hankel functions. The metric 
has the 10 symmetries of deSitter spacetime, and a solution having these symmetries is 

f n = (2Va(t) 3 )- 1 ^ g Hl{v). (7) 

with 


v = kH 1 exp (—Ht) (8) 

i/ = (9/4 - m 2 /H 2 ) 1/2 (9) 


Note that v decreases with increasing t. Alternatively, one could replace Hl(v) in (J7J) by H 2 (it) 
with u = —v , where the variable u increases with t. (See Q] , Sec.2.10, for a discussion of deSitter 
spacetime and the Bunch-Davies vacuum state.) 

As already mentioned in Section Q] the linearized Einstein equations for a gravitational wave 
field in the Lifshitz gauge reduce to a pair of massless minimally-coupled scalar field equations, 
one equation for each polarization of the gravitational wave. It is now clear that quanta of the 
gravitational wave field, i.e., gravitons, would be created by the expansion of the universe, in 
the same way that massless, minimally-coupled scalar particles would be created. For detailed 
quantization of linearized gravitational waves, see [10]. 

Now let us continue with our discussion of the minimally-coupled scalar field to summarize the 
method that I used to arrive at the number density of particles created from the initial Minkowski 
vacuum state (i.e., spontaneous particle creation) and from a non-vacuum state (i.e., stimulated 
particle creation). Since the reaction back of the created particles would be very significant in 
the early expansion of the universe, I chose to deal with an arbitrary expansion of the universe 
that began in a Minkowski spacetime and ended in a Minkowski spacetime (to within a constant 
rescaling of the coordinates) and that had continuous functions a(t), a(t), and a(t). 

From the theory of ordinary differential equations there are two linearly independent solutions 
of (J6l) . One can prove that the particle number at late times is an adiabatic invariant and that 
the solutions that approach positive and negative frequency solutions in the late time Minkowski 
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spacetime approach the asymptotic forms of the Liouville adiabatic solutions of the harmonic 
oscillator equation with slowly changing frequency. (This is also the form of the JWKB solutions 
for a slowly changing quantum mechanical scattering potential in one spatial dimension, with the 
particle energy above that of the scattering potential - but in the present case, the spatial coordinate 
x is replaced by the time coordinate t.) 

The adiabatic approximation of Liouville for the harmonic oscillator with a slowly changing 
frequency gives the following approximation to the solution of ([6]): 

hk(t) ~ (w fc (f)) -1/2 exp(±i J u k (t')dt'), (10) 

where 

cu fc (t) = \] ( k/a(t )) 2 + m 2 (11) 

This approximation is good in the limit that all time-derivatives of a(t) smoothly approach 0. 
The two solutions in m are linearly independent. Then the general solution of the second- 
order ordinary differential equation (|6|) can be written in the adiabatic approximation as a linear 
combination 

h k {t) ~ a k (w fc (t)) _1/2 exp(-i J Lo k (t')dt') 

+/3 k (u k (t))~ 1/2 exp(+i J u k (t')dt'), (12) 

where a k and (3 k are complex constants that must satisfy 

\a k \ 2 -\p k \ 2 = l ( 13 ) 

because of the conserved Wronskian or scalar product. We will assume that a(t) and all its time- 
derivatives are smooth and well-defined. 

In order to determine the expectation value of the number of particles that are created as 
a result of a smooth expansion of the universe, we take a(t) such that lim^-,*, a(t) = a±, and 
lirn^oc a(t) = a 2 , where a(t) and all its time-derivatives exist and are well-defined and continuous, 
and 0 < ai < 02 - 

We can take the solution of ([HI) in the early-time Minkowski spacetime such that |afc| 2 = 1 and 
|/3fc| 2 = 0. Then the expression for the quantized field in ([4]) is the standard expression for the 
quantized field in the early time Minkowski spacetime, for which a(t ) = a\. 

In the late-time Minkowski spacetime, when one uses the late-time asymptotic form in (1121k 
and separates the positive frequency and negative frequency terms in the late-time expression for 
quantized field (f> in (|^|). one finds that the field (f> at late times can be rewritten in the form: 
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<£ = S { *) + B \a*k (*» *)} 


(14) 


where 

<?£ = (21/a 2 3 )- 1 / 2 e^exp(-i J* uj k {t') dt'). (15) 

Here k l = 27rn*/L with n* an integer, and k = |fc|, and 

= + Pk —k' 

This transformation corresponds to the creation of particle-antiparticle pairs of total momentum 
0 from the vacuum, as required by conservation of momentum. This linear transformation is an 
example of a Bogoliubov transformation, as I explained in Q] on page 1061, in connection with 
particle creation by expanding universes. Later, analogous Bogoliubov transformations involving 
sets of creation and annihilation operators were used by Pulling [15] to describe particles in the 
spacetime (Rindler spacetime) corresponding to a set of accelerating observers and by Hawking 
161 ] to describe the distribution of particles created by a black hole. 

In the expanding universe under consideration here, the expectation value of the number of 
particles created by the expanding universe from the early-time vacuum state |0), present at late 
times is 


(0|Ht S | 0 ) = \/3 k f 


(16) 


(Recall that we are using the Heisenberg picture to describe the dynamics, so the state vector 
does not change with time.) Hence, if (3 k is non-zero, there is particle creation as a result of the 
expansion of the universe between the initial and final Minkowski spacetimes. 

Except in very special cases, quanta of the minimally-coupled scalar held are created by the 
expansion of the universe. From our earlier discussion of the quantized linearized-Einstein grav¬ 
itational held in the Lifshitz gauge, the same should be true for the creation of gravitons by the 
expansion of the universe. For the isotropically expanding FLRW universes, the particles are cre¬ 
ated in pairs to conserve the 3-momentum. I also calculated the detailed probability distribution 
of the created particles as a function of momentum. I used similar methods in Q to obtain the 
probability distribution of particles created by a Schwarzschild black hole. 

I also considered the creation of particles if the initial state is not the Minkowski vacuum state, 
but instead has a non-zero density of minimally-coupled scalar particles present at early times. 
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More generally, I considered a mixed state described by a statistical density matrix operator p stat 
in the initial Minkowski spacetime. I showed that there would be stimulated creation of these 
bosonic particles by the expansion of the universe. The details of the calculation are described in 
4], and in j(| on pp. 1063-1064. 

Taking an isotropic distribution of particles initially, and an FLRW expansion, the average 
number of particles present in mode k at late times in the final Minkowski spacetime is 


{N^)iate — Tr[p s t a tB^B^\ — (N~) ear i y + |/3fc| 2 (l + 2(N~) ear i y ) 


(17) 


The stimulated creation of bosons is evident in the factor involving |/3*.| 2 (1 + 2 (Nr) ear i y ). The 
average number of these spin-0 bosons present at late times is increased by stimulated creation 
resulting from the presence of a non-zero average number of these particles at early times. I also 
showed in (4], and in [7J] that the opposite is true in the case of fermionic particles created by the 
expanding universe. 

But as we shall see in Section m, when we discuss the rapid rate of increase of \/3k\ 2 during 
the exponential stage of the expanding universe, the initial number of particles present in the 
coordinate volume L 3 can greatly influence the creation rate and the magnitude of the number of 
particles present at late times. 

A further treatment of the stimulated creation of scalar gravitational perturbations during 
inflation is given by Ivan Agullo and me in |l9, 20j] where we discussed its possible effect on the 


CMB bispectrum. In 


211 ] with J. Navarro-Salas, we considered the CMB trispectrum. 


III. EXPONENTIALLY EXPANDING FLRW UNIVERSE 

Next we show how the above equations, (11611 and (117ft . for the expectation value of the number 
of minimally-coupled scalar particles created in each mode k determine the time-dependence of 
|/3fc| 2 in an exponentially-expanding, spatially-flat, FLRW universe. 

Let us start with a quantized minimally-coupled scalar field in an initial flat spacetime invariant 
under the 10-dimensional Poincare group. Suppose that the quantized scalar field is in its vacuum 
state having the 10 Poincare symmetries. Let us assume that our universe undergoes a spontaneous 
fluctuation (with a(t), a(t), and a(t) continuous) that eventually takes it to an exponentially- 
expanding, spatially-flat universe with a(t) oc exp (Ht), with H constant. After a time, suppose 
the quantized scalar held in this exponentially expanding deSitter universe has settled into a state 
having the 10 deSitter symmetries. Because we are working in the Heisenberg picture, the state 
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vector does not change with time, so it is the state |0) that is annihilated by the operator Ar, 
the annihilation operator for the scalar quanta in the original early-tinre Minkowski spacetime. 
However, the state |0) is not annihilated by the operators B g that annihilate particles in the late¬ 
time Minkowski spacetime. Recall that In the Heisenberg picture, the quantized field operator cj) 
evolves with time, and the B^ in the late-time Minkowski spacetime are obtained from the held 
operator at late-times. The particles that are present in the late-time Minkowski spacetime are the 
real (as opposed to virtual) particles that were created during the exponential de-Sitter portion of 
the expansion and thus survived the gradual smooth slow-down of the expansion to approach the 
late-time Minkowski spacetime. 

We will show from this and the equation for particle creation that |/3^(t)| 2 must be proportional 
to exp (Ht)/ exp(i7fi), where t\ is approximately the time that deSitter symmetry was reached by 
the particles being created, and t > t\. We have from (flfil) . in the Heisenberg picture, for particle 
creation by the expanding universe from the Minkowski vacuum state in the initial Minkowski 
space that 

(0|Rt^| 0 ) = 

Suppose that the expansion a(t) approaches deSitter, with a(t) oc exp (Ht) smoothly, with at least 
a(t), a(t), and a(t) continuous. (This is so that the minimally-coupled scalar field equation is 
well-defined at all times.) Here k l = 2'Kri 1 /L with ri 1 an integer, and k = \k\, and 

B%{h) = a k (ti) Aj: + 

This transformation corresponds to the creation of pairs of total momentum 0 from the vacuum, 
as required by conservation of momentum. We also have as a consequence of the field equation 
that, as in (fl3l) . 

|ct!fe(ti)| 2 — |/3fc(H)| 2 = 1, 

which is the condition for the commutation relation of boson creation and annihilation operators 
to obey Bose-Einstein statistics, i.e., to satisfy the correct commutation relations. 

After deSitter invariance has been reached, we expect that the number of particles present per 
unit physical volume ( La(t ) 3 ) is constant for t > t\. Hence, if t > t\, it follows that 


\Pkit)\ 2 / a{tf 
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is constant. Therefore, during the period for which inflationary equilibrium has been reached, we 
must have 


\/3k(t2)\ 2 /a(t 2 ) 3 = \Pk(ti)\ 2 /a(ti 


which implies that 


\Pk(t)\ 2 oc (exp {Ht)) 3 , 


(18) 


L 


after inflationary equilibrium has been reached. 

This result in (1181) was confirmed numerically by Glenz and Parker m by means of an explicit 
numerical integration, in which we took an expansion scale factor a(t) that started expanding 
from flat spacetime, eventually becoming an exponentially expanding inflationary universe that 
expanded for more than 50 e-foldings and then gradually approached a flat spacetime at late times. 
We joined the three segments of the expansion such that a(t), a(t), and a(t) were continuous at 
the two joining points. We used the known analytic solutions of the quantum field differential 
equations in each of the three segments of the expansion of the universe. As a check, we confirmed 
that (fl3l) was satisfied to at least 500 significant figures at all times throughout the expansion of 
the universe. 

As an interest ing aside, let me mention that the form of the parametrized scale factors a(t) 


that we used in 


18] to smoothly join the exponential inflationary expansion of the universe to 


Minkowski spaces at early and at late times are instances of the same four-parameter scale factor 
that I had used in a 1976 paper in Nature[17] to demonstrate that an expanding universe can 
create free particles from vacuum, with the created particles having a thermal distribution in the 
absence of any thermalizing interactions. 

The time-dependence of \/3k(t)\ 2 in (|18D also follows in the case when an initial distribution of 
particles is present, as in (fT7|) . This is clear because the average number of created particles that is 
proportional to the average number, (Nr) ear i y , of particles already present at early times in (fI71) is 
(1 + 2\/3k\ : ){N^) ear i y ~ 2\f3k\ 2 {N^) ear i y , so the argument given above leads in this case also to the 
same time-dependence, (fT8j) . of |/3fc| 2 . However, we see that the average number of particles created 
during the exponential expansion can be greatly increased by the stimulated creation of particles, 
if there were a significant average number (N^) ear i y of these bosonic particles already present in 
the physical expanding volume when the exponential stage of the expansion first began. 
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